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The detection of a nuclear spin in an individual molecule represents a key challenge in physics 
and biology whose solution has been pursued for many years. The small magnetic moment of a 
single nucleus and the unavoidable environmental noise present the key obstacles for its realization. 
Here, we theoretically demonstrate that a single nitrogen-vacancy (NV) center in diamond can be 
used to construct a nano-scale single molecule spectrometer that is capable of detecting the position 
and spin state of a single nucleus and can determine the distance and alignment of a nuclear or 
electron spin pair. In combination with organic spin labels, this device will find applications in 
single molecule spectroscopy in chemistry and biology, such as in determining protein structure or 
monitoring macromolecular motions and can thus provide a tool to help unravelling the microscopic 
mechanisms underlying bio-molecular function. 



Introduction. — As single nuclei have exceedingly small 
magnetic moments, a large ensemble (typically 10 18 nu- 
clear spins) is necessary to obtain an observable signal 
employing methods such as magnetic resonance spec- 
troscopy (NMR). As a consequence, chemical and bi- 
ological processes have thus usually been tracked with 
ensemble measurements, which only provide ensemble 
averages and distribution information. Single-molecule 
studies can instead allow one to learn structural infor- 
mation and time trajectories of individual molecules free 
of natural disorder PHI]- The detection of a single nu- 
cleus can thus provide various new possibilities for single 
molecule spectroscopy. Furthermore, single nuclear spins 
have long coherence time due to their weak coupling with 
the environment which makes them promising candidates 
for a qubit or a quantum register. For example, it has 
been proposed to engineer nitrogen nuclear spin in the 
molecule of 14 N@C6o [3 E] or phosphorus donor [7] as 
a qubit candidate. For these purposes efficient readout 
of the nitrogen nuclear spin state of a single molecule is 
crucial. 

Nitrogen-vacancy (NV) defects in diamond have been 
used to construct ultrasensitive nano-scale magnetome- 
ters [8 13], and have led to interesting applications in 
nano imaging [14] [15] and biology as well [T6HT8] . N V 
centers in diamond benefit from long coherence times 
at room temperature and highly developed techniques 
for coherent control and optical readout of their elec- 
tron spins. A single NV centre can be used to detect 
single external electron spin [19] and a proximal nuclear 
spin (within a few atom shells) [20-22 . It has also been 
shown that it is possible to detect a strongly coupled nu- 
clear spin pair [23] and measure noise spectra [24] [25] 
by applying dynamical decoupling pulses. The scheme 
nevertheless is not applicable for sensing single nuclei 
[23H25] , as the effect of a single nucleus is completely 
removed by the decoupling pulses and is featureless in 
the dynamics of NV spin. As a consequence the detec- 
tion of individual nuclear spins in the presence of realistic 



environments remains an unsolved challenge. Key obsta- 
cles originate firstly from the requirement that coherence 
times are sufficiently long in order to observe the effect 
of single nuclei on NV spin. The second key obstacle lies 
in the difficulty of distinguishing a distant nuclear spin 
from other environmental nuclei that couple to the NV 
spin. 

In this work, we address both challenges by continu- 
ously driving a single NV spin and use it as a probe to 
measure a specific transition frequency of the target sys- 
tem. The role of continuous driving is two-fold: firstly, 
it decouples the NV spin from the unwanted influence of 
a spin bath (which is particularly strong for NV centers 
located close to the surface [19]) to achieve sufficiently 
long coherence times [26-28 ; secondly, by changing the 
Rabi frequency of the external driving field we can tune 
the NV center to match the target frequency and thus 
selectively enhance the sensitivity for this specific fre- 
quency and thus single out the target system. With this 
scheme, we are able to determine the position of a sin- 
gle nucleus in the presence of a realistic environment. 
We demonstrate how to use this mechanism to imple- 
ment quantum non-demolition(QND) measurement of a 
single nitrogen nuclear spin state in a cage molecule of 
fullerene ( 14 N@Ceo)- This is achieved by exploiting the 
fact that the flip-flop between the NV spin and the target 
system will either be allowed or prohibited dependent on 
the nuclear spin state for a specific NV spin initial state. 
Our proposed detector can be applied straightforwardly 
to the determination of the distance and alignment of a 
spin pair. We show that the present model can find appli- 
cations in single molecule spectroscopy with organic spin 
labels. It can also be used to witness the creation and 
recombination of charge separate state in radical pair re- 
actions. We expect that the present diamond-based sin- 
gle molecule probe will have more potential applications 
in determining protein structure and monitoring chemi- 
cal (biological) processes, see the examples of molecular 
motors, ATPase and RNA folding in [SJHSI]. 
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FIG. 1. (Color online) Model of a NV center spectrometer. 
We apply a continuous field to drive the electronic transition 
\m s = 0) <H- \m s = — 1). When the Rabi frequency Q is on 
resonance with one specific transition frequency of the target 
system, the flip-flop process will happen and lead to the decay 
of the NV dressed state population. 



Basic model of NV center spectrometer — The probe 
in our model is a nitrogen-vacancy (NV) center spin in 
diamond, whose ground state is a spin-1 with a zero- 
field splitting of 2.87GHz. By applying an addition mag- 
netic field, one can lift the degeneracy of \m s = —1) and 
\m s — +1) and thus allow selective driving with contin- 
uous microwave field of one specific electronic transition, 
e.g. \m s = 0) \m s = +1) (henceforth denoted by |0) 
and |+1) respectively) . Within the |0), |+1) subspace 
the Hamiltonian of the driving field can be written as 
Hnv — Qo~ x , where a x is spin-^ operator and the eigen- 
states represent the dressed states of the system. We re- 
mark that different driving schemes may be used depend- 
ing on the properties of target systems as we will discuss 
later. Our goal is to use such a NV dressed spin as a 
probe to detect a specific frequency in the target system. 
The interaction between the NV spin and another spin 
is H NV _ S = ^2n9n [3 (S ■ f) (I N ■ r) - S ■ I N ], where S, 
I/v are the NV spin and target spin operators, = 
-(hfjLoj N j e )/(A7rR%), 7e and j N are the gyromagnetic 
ratio of electron spin and target spin respectively, f is 
the unit vector that connects NV center and the target 
spin. The large zero-field splitting leads to energy mis- 
match which prohibits direct NV spin flip-flop dynamics 
and allows for the secular approximation to simplify the 
NV-target spin interaction to 

H N v-s = S z J2dN [3r z (r x I x N + r y l y N ) + (3r 2 z - l) 1^] 

N 

(i) 

The coupling operator S z leads to dephasing type inter- 
action in the bare basis but induces flips of the NV spin 
in the dressed states basis. The flip-flop process will be 
most efficient when the Rabi frequency i.e. the energy 
splitting of the dressed states, matches the the transi- 
tion frequency of the system spin that we wish to probe. 
Therefore, one can initialize the NV spin in one of the 
dressed states and tune the Rabi frequency to determine 
the transition frequency of the target system by mea- 



suring the probability that the NV spin remains in the 
initial state. The target system distance can be inferred 
from the flip-flop rate of the NV center if the magnetic 
moment is known. 

Measure position of a single nucleus. — We first ap- 
ply our idea to detect the position of a single nucleus. 
This is particularly interesting in biology, for example 
in determining where certain important nuclei (e.g. 14 N, 
31 P) are located in the protein complex. We consider the 
combined system of NV spin under driving and a target 
nuclear spin with the following effective Hamiltonian 



- u N l • h + ^ 



(g^rj+l) a z (i • h) , (2) 



where o z is the spin-^ operator, and the hyperfine vec- 
tor h is determined by the position vector as h x = 
3r x r z / v / 3r 2 + 1, h y = 3r y r z /y/3r 2 z + 1 and h z = (3r 2 z - 
1)/ \/3r 2 + 1. The effective nuclear spin Larmor fre- 
quency ujn is determined by ujn = jNB e with B e = 
B(0,0) - \ (g\/3r 2 z + l) h. We choose the Rabi fre- 
quency Q on resonance with the nuclear Larmor fre- 
quency to match the Hartmann-Hahn double resonance 
condition as Q = [32]. After preparing the NV 
spin in the initial state |+) = ^ (|0) + |+1)) and as- 
suming that the nuclear spin is in the thermal state 
namely /?at(0) = \l at room temperature, the signal 
S(t) = (+|p(t) |+) at time t is 



S(t) 



1 + cos (Jt) 



(3) 



where J = \ {^g^3r 2 z + 1^ [x — h-b^ is the effective spin 
flip-flop rate. The effect of the 13 C spin bath in diamond 
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FIG. 2. (Color online) Measure the position of a single nu- 
cleus 31 P in a 1 H3 31 P04 molecule at a distance of d — 5nm 
from the NV center, (a) The signal S(t) measured at time 
t — 3ms for different magnetic field directions (0, <f>) with the 
application of radio frequency driving (20kHz) on resonance 
with the Larmor frequency of ^^H. The driving amplitude on 
NV spin is tuned to be on resonance with the Larmor fre- 
quency of 31 P at uj — j 31 pB = 500kHz. The hyperfine vector 
is in the direction of O = 68.233° and O = 93.841°. (b) The 
signal S as a function of time t. The magnetic field direction 
is chosen to be orthogonal to the hyperfine vector. The esti- 
mated value of J is 0.2041kHz as compared to the exact value 
0.2065kHz. 
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can be greatly suppressed if the 13 C environmental spec- 
tral density is low at the Rabi frequency Q of the driving 
field. We can use an ultra pure diamond sample, and thus 
the intra-interaction between 13 C spins is weak. The rel- 
evant frequency of the spin-bath is mainly determined 
by the applied magnetic field, namely uj = ji3 C B. If 
the mismatch between Q and uj is sufficiently large, the 
NV spin will be effectively decoupled from the 13 C spin 
bath. There is a possibility to grow isotropically pure 
12 C diamond [33 . The sensitivity for measuring the tar- 
get nucleus is then determined by the relaxation time T\ 
which is very long for an NV spin (T\ ~ 10 ms at room 
temperature) [34] . 

We demonstrate our idea by showing how to use a 
NV spin to measure the position of 31 P in a single 
molecule 1 H3 31 P04. Due to its half-integer nuclear spin 
and high abundance of 31 P, NMR spectroscopy based 
on phosphorus-31 has become a useful tool in studies 
of biological samples. Phosphorus is commonly found 
in organic compounds, coordination complexes and pro- 
teins, such as phosphatidyl choline which is the major 
component of lecithin. In 1 H3 31 P04, 31 P interacts with 
a few adjacent hydrogen atoms, which captures generic 
features of the protein. Thus it can serve as a paradig- 
matic example for practical applications of our model. 
We apply a magnetic field such that the Larmor frequen- 
cies of 31 P, X H and 13 C are different. The Rabi frequency 
of the NV spin is tuned to match the (effective) Larmor 
frequency of 31 P and is sufficiently strong to eliminate the 
effect from the 13 C spin bath, see SI for details. The line 
broadening caused by the adjacent X H is relatively large 
when compared with the coupling between the NV spin 
and 31 P. To overcome this problem, we propose to con- 
tinuously drive 1 H nuclei, which effectively averages out 
the line broadening. The difference of Larmor frequen- 
cies among different nuclei allows us to selectively drive 
1 H while not affecting the other nuclei. With this, we can 
obtain information of the hyperfine vector (and thereby 
the relative direction of 31 P), see Figga). One can then 
measure the flip-flop rate, from which the distance from 
31 P to NV center can be obtained from Eq.([3|, see SI 
for more details. We remark that shallow (below 5 nm) 
implantation was realized experimentally [35] . 

Quantum non- demolition measurement of a single nu- 
clear spin state. — In the above model, the flip-flop pro- 
cess happens under the resonance condition but also re- 
quires that the nuclear spin state is complementary to NV 
spin state. As nuclear spins have long coherence time and 
can serve as robust qubit we would like to show that it is 
possible to perform QND measurement on the system of 
14 N@C6o i-e. a nitrogen atom in a Ceo cage. The 

14 N@C6o molecule has an electron spin-| coupled to the 
14 N nuclear spin-1. The hyperfine interaction is isotropic 
and the spin Hamiltonian is given by 
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FIG. 3. (Color online) (a) Scheme for quantum non- 
demolition measurement of the nitrogen nuclear spin state 
in 14 N@C6o with a NV center at a distance of 8nm. For the 
nuclear spin state |mj = 0), the allowed electron transition 
frequency in 14 N@C6o is uj e and will not be on resonant with 
uJnv, while for \mi = +1), the resonant condition is satis- 
fied when uj e + a — (2Q 2 +cje) 1 ^ 2 - The signal S and the 
nuclear spin state fidelity as a function of time t for the nu- 
clear spin state \mi = +1) (b) and \mi = 0) (c). The Rabi 
frequency Q, is chosen as Q = 70.7MHz corresponding to the 
resonant frequency Qnv — 316.2MHz for the nuclear spin 
state | mi — +1). The applied magnetic field is 300MHz. The 
electron spin of 14 N@Ceo is not polarized, and its initial state 
is approximated by the maximally mixed state. 



where uj e = — 7 e 5, wn = JnB and the quadrupole 
splitting is Aq = 5.1MHz, the hyperfine coupling is 
a = 15.88MHz. We can encode a qubit in the nuclear 
spin state |0) 7 = |raj = 0) and |1) 7 = \mj = +1). As 
the electron and nuclear Zeeman splitting exhibit a large 
mismatch due to a strong magnetic field, the main effect 
of hyperfine coupling is to modulate the electronic tran- 
sition frequencies dependent on the nuclear spin state. 
When the nuclear spin state is \mj = 0), the allowed elec- 
tronic transition frequency is cj e , while for the nuclear 
spin state \mi = +1), it is uj e -\-a. We note that under an 
additional magnetic field, the energy separation of NV 
spin state \m s = +1) and \m s = —1) is A = 2cj e , which 
means that we will inevitably drive both electron spin 
transitions (namely |0) O |=bl)) if the driving amplitude 
£1 is around u e in order to satisfy the resonance condi- 
tion. We thus apply a driving field at the frequency ujo 
with the detuning ±A for two allowed electronic transi- 
tions. The effective electronic transition frequencies be- 

2 i ,.,2A 1 /2 o , _ _ / oo2 , / .,2\ 1 /2 



H s = oj c S z n + uj n I z n + Aq (1 



N) 



aS 



N • J- TV 



(4) 



come uj\ = (2fl 
SI for details. 



0J 2 e ) 1/2 and uj 2 = 2(2Q 2 



, see 
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We prepare NV spin in the dressed state \D) = 
(SI |+1) + u e |0) - n |-1)) /(2Q 2 + a;2)i/2 and then tune 
the Rabi frequency Q to be on resonance with the allowed 
electron transition frequency in the system 14 N@C6o cor- 
responding to the nuclear spin state |raj = +1), namely 

1/2 

(2VL 2 + cjg) = uo e + a. It is important to note that we 
do not require the electron spin of 14 N@C6o to be po- 
larized and assume it is in the maximally mixed state, 
i.e. p e = | at room temperature. After a time t, if 
the nuclear spin is in the state |raj = 0), the resonance 
condition will not be satisfied and the NV spin will re- 
main in its initial state \D); if the nuclear spin is in state 
| raj = +1) the population of state \D) will change with 
time. The nuclear spin state remains unaffected through- 
out. Thus the readout represents a non-demolition mea- 
surement that can be repeated. The feasibility of this 
idea is verified by our numerical simulation, see Figj3] 
It can be seen that the required time for one readout is 
much shorter than the coherence time T2 of the electron 
spin in 14 N@Ceo (which is 20/is at room temperature 
[5]) when the distance between 14 N@Ceo and NV center 
is 8nm. We remark that the present mechanism can be 
exploited to polarize nuclear spins, see SI for details. 
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FIG. 4. (Color online) (a) Resonance of a single 1 H20 
molecule with an additional magnetic field Ji^B = 500kHz in 
the x directions. The signal S is measured at time t = 0.6ms. 
(b) Resonance of a pair of nitroxide spin labels with a distance 
of 8nm. The magnetic field along the vector which connects 
two spin labels. Continuous driving field with Rabi frequency 
of 20MHz is applied on the spin labels. 

Measure distance and alignment of a spin pair. — For 
the tracking of molecular conformations or the determi- 
nation of molecular structure it can be of advantage to 
determine the distance between two specific spins. To 
this end, we consider a target system which consists of 
two spins interacting with each other. To determine their 
intra-distance R and the orientation of the alignment vec- 
tor f, we apply a strong magnetic field ujn = 7JV-B, as 
the energy shift will depend on the relative orientation of 
distance vector and applied magnetic field. In the mean 
time make the NV spectrometer will be operated with a 
strong driving field to decouple from the 13 C spin bath in 
diamond. We tune the amplitude of driving on the NV 
spin around the Zeeman energy a; at, so that the dominant 



resonant frequencies are O = ujn± 



l-3(f-b(0,<P)f 



see SI for more details. The exact values of g and f can 



be inferred by applying magnetic fields in nine directions 
and measuring the resonant frequencies respectively (see 
SI for more details). 

As an example, we apply our idea to measure the dis- 
tance between two hydrogens in a water molecule. In 
Fig(4j we plot the resonance frequencies with a magnetic 
field applied in the x direction. The results for other 
magnetic field directions can be found in supplementary 
information. The 1 H20 molecule is 5nm from the NV 
center, and the vector that connects two hydrogen atoms 
is described by (r,6> o ,0o) = (0.1515nm, 118.2°, 288.85°). 
With the resonance frequencies, we can infer that g = 
34.684kHz, which implies an the intra-molecular dis- 
tance of d = 0.1518nm and the alignment direction is 
((9,0) = (118.29°, 288.82°), see SI for more details. The 
difference from the exact values is due to the slight Lar- 
mor frequency shift induced by the NV spin. If the spin 
pair is rotating randomly for example in solution, the 
resonance frequency will only dependent on the intra- 
distance. 

The above idea can be combined with spin labels and 
have potential applications in chemistry and biology. 
Spin labels are organic molecules with a stable unpaired 
electron [36]. They can be attached to the protein (co- 
valent or as a ligand) via a functional group. It was pro- 
posed to use a spin label to amplify the sensitivity of a 
diamond manometer [12]. Electron spin resonance (ESR) 
on an ensemble based on spin labels has widely been used 
as a spectroscopic ruler to determine protein structure 
and monitor macromolecular assembly processes. How- 
ever, it is very hard to go beyond a distance of 5nm be- 
tween spin labels, due to e.g. inhomogeneous line broad- 
ening. We consider the widely used nitroxide spin labels 
and show that the resonance linewidth (which mainly 
depends on the coupling strength between NV spin and 
spin labels) assisted by dynamical nuclear polarization 
and continuous drivings can be smaller than the reso- 
nance frequency splitting for a pair of spin labels with 
a distances larger than 5nm (more details are present in 
supplementary information), see Fig. 4 (b) for an example 
of 8nm. In a similar way, one can monitor the creation 
and recombination of radical pairs, e.g. in reaction cen- 
ters, see supplementary information for more details. 

Summary. — We have proposed a scheme to construct 
a nano-scale single molecule spectrometer based on NV 
centers in diamond under continuous driving. This spec- 
trometer is tunable by changing the Rabi frequency. We 
demonstrate its application to the detection of a single 
nucleus, including its position and spin state. The idea 
can also be used to measure the distance and alignment 
of a spin pair. This opens a novel route to determine the 
structure of proteins and monitor conformational changes 
and other processes of relevance to biology, chemistry and 
medicine. We expect that our result and its extension can 
greatly enrich the diamond-based quantum technologies 
and their applications in chemistry and biology. 
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SUPPLEMENTARY MATERIAL 
Continuous dynamical decoupling of the 13 C spin bath 

The NV center spin is coupled with a 13 C spin bath, the interaction is described as follows 

Bsv-»C = ^ E 1? P • I™ " 3 (S • r) (I m • f )] , (5) 

where R m is the distance from NV center to the nuclei, and f is the unit vector that connects NV center and nuclei. 
The Hamiltonian of the spin bath itself is 

tfisc = - E ^ B • J » + E ^ ^ ■ - 3 ^ • f ) • ^ ■ ^ 

n n>m " bn 

where R mn is the distance between two nuclei, and r is the unit vector that connects two nuclei. The mass density of 
diamond is p = 3.5g • cm -3 , so the number density of carbon is n = 3.5/(12 x 1.66 x 10 -24 )cm -3 = 176nm -3 . For a 
sample with 0.01% 13 C, the number of 13 C in a sphere with the radius 4nm centered at NV center is 176 x ^ x 4 3 x 
0.01% « 5. As we want to single out the effect of the target system, we apply continuous driving filed on NV spin. 
Here, we use exact numerics to show that our continuous dynamical decoupling is very efficient for the parameters 
that we use for the examples in the main text, see Fig JH] We consider a diamond sample with 8 13 C spins in a 4 
nm sphere. We remark that, in a similar way, continuous dynamical decoupling will also suppress the noise from the 
surface of diamond if it is due to spins that are different from the target spin. 
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FIG. 5. (Color online) Decouple NV spin from the 13 C spin bath with continuous dynamical decoupling. The signal S(t) = 
(+|p(£) |+) of NV spin as a function of time t. (a) In the example of 1 H| 1 P04, the magnetic field is j 31 pB = 500kHz 
(i.e. Ji 3 qB = 310.6kHz), and the driving amplitude is Q — 500kHz. (b) In the example of 1 H20, the magnetic field is 
jillB = 500kHz (i.e. Ji 3 qB = 125.7kHz), and the driving amplitude is Q = 400kHz. 



Measure a single 31 P position in the molecule of 1 H| 1 P04 



The effective hyperfine interaction between NV spin and the target spin is 

H NV _ S = gSz [(3r* -1)I Z N + 3r z (r x F N - r y \ y N )} . 



(7) 



We introduce the hyperfine vector /i(0o,</>o) 7 which is determined by the vector r that connects the NV center and 
the target spin as 

3r?-l 



cos (0 O ) 
sili (0 O ) cos 0o = 
sin (0o ) sin 0o = 



3r x r z 
3r y r z 



The effective Hamiltonian is thus rewritten as 

H E = nS x - ( lN B) I • b (0, 0) + (gy/Zr* + l) S z 



l-ft(flo,0o) 



We note that S z = |+1) = \l + o z where o z is the spin-^ operator, thus the above Hamiltonian is 



H E = nS x - lN I • B e (0 e , <t> e ) + X - (gy/3r* + l) a z 



l'h(0 ,<fo) 
B 



(8) 
(9) 
(10) 

(11) 

(12) 
We 



where the effective nuclear spin Larmor frequency is determined by B e = B(0,</>) — \ {^g^f^r\ + 1^ /z(0q,0o)- 
choose the driving amplitude ^o on resonance with the Larmor frequency of the target nuclear spin namely Q = 7Ar£> e . 
To measure the position of the target spin, we first prepare NV spin in the initial state |+) = ^= (|0) + |1)) and assume 

that the nuclear spin is in the thermal state which is well approximated by pjv(0) — If the Zeeman energy is much 
larger than the coupling (a; at = O ^> q\ the flip-flip between the NV spin and the nuclear spin is very unlikely, and 
the flip-flop is the dominant process. After time t, we measure the probability that the NV spin remains in the state 
|+) which is a function of the coupling strength J as follows 



S(t) 



1 1 + cos (Jt) 

2 + 4 ' 



(13) 



From the signal S(t) one can infer the value of J that depends on the distance and direction of the nuclear spin to 
NV center 



J 



i 



i-h(0 o ,4> o )-b(0 e ,4> e ) 



(14) 
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The spin flip-flop rate J as in Eq.(14) depends on the direction of the effective magnetic field (0 e ,0 e ). In the case 
that the applied magnetic field is much larger than the nuclear energy shift induced by NV spin, the direction (0 e , <fr e ) 
can be approximated well by the applied magnetic field direction 6(0,0). It can be seen that if and only if the 
magnetic field direction (0, (j>) is in parallel with the hyperfine vector, namely = 0o and <p = 0o, the effective flip-flop 
rate vanishes J = even though the Hartmann-Hahn resonance condition is still satisfied. If the magnetic field is 
orthogonal to the hyperfine vector, the flip-flop rate is maximal as J m = \ ^y / 3r| + 1^ . By choosing a series of 

magnetic field directions b(0, 0), one can obtain both the values of the coupling strength g and the direction h(0Q, 0o), 
from which one can obtain the direction of the target single spin (r xi r yi r z ) with respect to the NV center as follows 



r 2 




2 + cos 2 (9 //2 + cos 2 <9 \ 1 - cos 2 O 



T. 



y 



6 V V 6 / 9 

tan0 o . (15) 



The relative sign among r Xl r yi r z can be determined from Eqs.(8-10). 

We remarked that the fluctuation in the driving field (CI) will limit the measurement precision. One way to overcome 
this problem is by using the dressed qubit from concatenated dynamical decoupling as propose in pQ (instead of a bare 
spin qubit) as a probe. For example, we can use a second-order sequential driving scheme with = Cli cos(ut)a x and 
Hd 2 = ^2 cos [(uj + Ct)t] (j y , and get a qubit encoded in the second-order dressed states. In the second order interaction 
picture, we can tune the Rabi frequencies to satisfy the resonant condition Ct\ + ^2 = ^at, and the coupling operator 
of NV spin with the target system becomes 

(Jz^\ exp \i(Cl x + CL 2 )} \+ y ) (- y \ + h.c. (16) 

The energy gap of the dressed qubit {1+^) , \ — y )} is more robust against the error in the driving fields, see [I] for 
more details, and can thus make our scheme work with a higher accuracy. 

For 31 P and 1 H, the coupling strength is g = 48. 6(6. 075)kHz for the distance R m = 0.1(0.2)nm, which is relatively 
strong as compared the coupling strength between NV spin and the target 31 P spin. The sole information about 
the position of 31 P is blurred by such a strong interaction, see Fig{6] (a). We can apply a strong magnetic field such 
that (71 jj — 7sip)5 ^> g m , which suppresses the flip-flip (flip-flop) process between the target nuclear spins and its 
neighbouring hydrogen nuclei. The effective Hamiltonian of the target system thus becomes 

H s = — 7atB • lN - 7m B • Im + E 9m (1 - 3r z 2 ) I z m I z N (17) 




FIG. 6. (Color online) Measure the position of a single nucleus 31 P in a molecule of 1 H3 31 P04 at a distance d = 5nm to the 
NV center, (a) The signal S(t) measured at time t = 3ms for different magnetic field direction as described by (0, cj>) without 
radio frequency driving on resonance with the Larmor frequency of ^^H. The driving amplitude on NV spin is Cl — 500kHz. (b) 
The signal S(t) measured at time t — 3ms for different magnetic field direction as described by (6, 0) with the application of 
radio frequency driving (the amplitude is 20kHz) on resonance with the Larmor frequency of 1 H. The driving amplitude on 
NV spin is tuned to be on resonance with the Larmor frequency of 31 P. The hyperfine vector is in the direction of #o = 68.233° 
and 0o 93.841°. 
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To overcome the line broadening (i.e. the last term in the above equation), we propose to drive the hydrogen nuclei 
continuously. Since the Larmor frequencies can be made sufficiently different for 31 P and X H, it is possible for us to 
selectively drive the hydrogen nuclear spins while not affecting the target 31 P nucleus due to a large detuning. This 
leads to the system Hamiltonian as follows 

H s = - 7 ivB • lN + • J2 1™ + E 9m (1 - 3r z 2 ) V m V N . (18) 

m m 

If the Rabi frequency is large enough, the interaction term 1^1^ can be effectively averaged out and eliminated. 
We remark that one can also drive the target 31 P nuclear spin which can also suppress the line broadening caused by 
hydrogen nuclei. 

In our numerical simulation, we apply a magnetic field such that the Larmor frequencies of 31 P, X H and 13 C are 
quite different, namely 500kHz, 1235kHz and 310.6kHz respectively. The driving amplitude of the NV spin is tuned to 
match the effective Larmor frequency of 31 P, thus its interaction with the 13 C spin bath is eliminated. The difference 
of Larmor frequencies also allow us to selectively drive the nuclei X H (with the driving amplitude 20kHz) while not 
affecting the other two types of nuclei. The distance between the nucleus 31 P and X H is about 0.2nm, so that the 
coupling constant is g w 6KHz. The line broadening is much larger than the coupling between NV spin and the target 
nucleus 31 P (which is smaller than 1kHz at a distance of 5nm). In Fig|6|a), it can be seen that if we do not drive 
the nuclei X H, the interaction between 31 P and X H will smear the measured signal and we can hardly get information 
about the position of 31 P. By applying radio frequency driving field on resonance with the Larmor frequency of X H, 
we can then clearly obtain information of the hyperfine vector (and thereby the relative direction) of 31 P with respect 
to NV spin, see Fig{6jb). 



Initial state preparation of a nuclear spin bath 



Nuclear spins have long coherence times and can be used as robust quantum register [2] . It has been shown how to 
initialize the state of a nuclear spin proximal to the NV spin [3] and prepare a spin environment by optical cooling 
[4]. In the following, we demonstrate that our idea can also be exploited to initialize a nuclear spin bath with the 
distance from the NV center as far as a few nanometers with high efficiency and fidelity. We assume that the initial 
state of the nuclear spin is Pat(0) = |. The procedure of one cycle consists of two steps: first one prepares the NV 
spin in the state |+) = A= (|0) + |1)) and then drives the NV spin for time r with the Rabi frequency on resonance 
with the Larmor frequency of the nuclear spin. We define the nuclear spin up and down state along the direction 
of the applied magnetic field. After one cycle the nuclear spin down state \\) population P± increases from Pj r (fc) to 
P±(k + 1) = P+(fc)+ 1+cos 2 (2Jr) [1 - Pi(k)]. The formula will become more complicated if there are other nuclei. In fact 
one can prepare the initial state of several nuclear spins with the same procedure. The time r is constrained by the 
spin-lattice relaxation time T\. To overcome such a limit, we can choose a small time r and then repeat the above cycle. 

The nuclear spin state |^) population P± after k cycles can be approximately written as P±(k) = 1 — \ 1 — (Jr) 




5 10 
t(ms) 



100 




5 10 
t(ms) 



FIG. 7. (Color online) The polarization M of 13 C nuclear spins as a function of time t = kr in a diamond sample with dilute 
13 C. (a) 6 13 C within 4 nm from the NV center, the Larmor frequency uj 13 q on resonance with the NV spin continuous driving 
amplitude Q as cji3 C = O = 50kHz. (b) 8 13 C with cji3 C = Q = 500kHz. The cycle time is r = 0.5ms. 
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We consider a sample diamond with dilute carbon-13 nuclei (e.g. 0.01% 13 C) and show how one can initialize the 
quantum state of a single 13 C, which may later be used for quantum memory. It can be seen from Fig JT] that the 13 C 
can efficiently be initialized to a very high polarization (M = P± — P^). As the 13 C is very dilute in such a diamond 
sample, the flip-flop process between different 13 C is very slow, and thus each 13 C is mainly directly polarized by the 
NV center. The required time is then dependent on the distance and alignment between the individual 13 C and the 
NV center. 



Quantum non-demolition measurement of the nitrogen nuclear spin state in N@C6o 

Our idea can be used to readout quantum state of a single nuclear spin. We have shown that it is possible to 
perform quantum non-demolition measurement on the system of 14 N@Ceo, which has a nitrogen atom in a Ceo cage. 
The molecule has an electron spin 5 = § coupled to the 14 N nuclear spin 1 = 1. The hyperfine interaction is isotropic 
and the spin Hamiltonian can be simplified as 

H s = uj e S z N + uj n I z n + A Q (I Z N ) 2 + aS N • 1^, (19) 

where uo e = — 7 e 5, ojn = JnB and the quadrupole splitting is Aq = 5.1MHz, and the hyperfine coupling strength 
is a = 15.88MHz. We here consider a qubit in the nuclear spin state \0) T = \mj = +1) and |1) 7 = |raj = 0). If the 
difference between the electron and nuclear Zeeman splitting uj e — ujn is large, the non-secular terms in the hyperfine 
coupling can be neglected, and thus the system Hamiltonian is written as 

H s = uj n In + oj e S N + aS z N l z N . (20) 

When the nuclear spin is in the state of |raj = 0), the allowed electron spin energies is equidistant with uo ei while for 
the nuclear spin state |raj = +1), it is uo e + a. By applying an additional magnetic field, the energy separation of the 
NV spin state \ms = +1) and \ms = — 1) is A = 2cj e , namely 

H NV = (c^o - u e ) |-1) <-l| + (uj q + u e ) |+1) (+1| . (21) 

This means that we will inevitably driving both spin transitions |0) <-> |±1) if the driving amplitude is around uo e in 
order to satisfy the resonant condition between NV spin and the electron spin in 14 N@C6o- Thus, we have to take all 
three sublevels of the NV spin into account. Our idea is to apply a continuous driving field at frequency ujo = 2.87GHz 

as 

H d = ncos(u t) (|+1) <0| + |-1) (0| + h.c.) . (22) 
In the interaction picture with respect to Hnv, we have 

Ht = (2tf + to 2 .) 1 / 2 (\u) (u\ - \d) (d\) , (23) 

where we have three dressed states as 

I") = fa+ I" 1 ) + 2? ?+ 1°) + l + ' ( 24 ) 

//_(_ -h z 

lD) = (2^+^2)1/2 (" l + + ^ 1°) " Q I" 1 )) ' ( 25 ) 

\d) = ^^(v 2 -\-l)-2V-\0) + \+l)), (26) 

where rj± = [(2Q, 2 + ujI) 1 ' 2 ± u e ] /fi. We note that the NV spin coupling operator S z with the other spins will induce 
transitions among the dressed states with two transition frequencies uj\ = + oj 2 ) 1 ^ 2 and UJ2 = 2(2^ 2 + uj 2 ) 1 ^ 2 . 
Therefore, we first prepare NV spin in the state \D) and then tune the Rabi frequency to be on resonance with 
the allowed electron transition frequency in the system 14 N@C6o corresponding to the nitrogen nuclear spin state 
| raj = +1), namely 

(2Q 2 +ou 2 e ) 1/2 =u e + a. (27) 

After time r, we measure the probability that NV spin remains in the state \D). If the nuclear spin is in the state 
| raj =0), it can be seen that the resonant condition would never be satisfied once we drive the NV spin (namely 
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FIG. 8. (Color online) Measurement of the nitrogen nuclear qubit state in 14 N@C6o with a NV center at a distance of 8nm. 
The signal S measured at time t = as function of the effective NV transition frequency ujnv = (2Q, 2 + oo 2 ) 1 ^ 2 for two 
nuclear spin state \mi — 0) (red) and |raj = +1) (blue). When \mi — 0), the allowed electron transition frequency in 14 N@C6o 
is uj e and will not be on resonant with ujnv- If \mi — +1), the resonant condition is satisfied when uo e + a = (2Q 2 + uo 2 ) 1 ^ 2 . 
The splitting of the resonant frequency stems from the virtual transitions between the electron and nuclear spin in 14 N@C6o 
caused by the non-secular hyperfine coupling terms. The applied magnetic field is 300MHz. The electron spin of 14 N@C6o is 
not polarized, and its initial state is approximated by the maximally mixed state. 



fi 7^ 0), NV spin will thus stays in the initial state \D)\ otherwise if the nuclear spin state is \mj = +1), the flip-flop 
process can happen and thus the state population of \D) will change with time at the resonant frequencies. In Fig(8j 
we see that there are actually three resonant frequencies. The splitting comes from the virtual electronic transition 
caused by the nuclei (~ a 2 /uo e = 0.84MHz). In the mean time, the nuclear spin state populations are not affected by 
the readout procedure, and one can thus implement repetitive measurement of the nitrogen nuclear spin state. In our 
numerical simulation, we apply a magnetic field \ j e B\ = 300MHz, and the corresponding Larmor frequencies of 13 C 
and 14 N are 114.75kHz, 32.97kHz respectively. The amplitude of driving on the NV spin corresponding to the central 
resonant frequency ujnv — 316.2MHz is Q — 70.7MHz, which is strong enough to suppress the effect of the 13 C spin 
bath. We remark that if we first polarize the electron spin to the state | — |) (e.g. by using NV center), it is possible 
to improve the readout efficiency. 



Measure distance between two hydrogen nuclei in a water molecule 

We consider a target system which consists of two spins interacting with each other. The system Hamiltonian is 
written as 

H s = -7ivB • (l Nl + I N2 ) + g [I Nl • Iiv 2 - 3 (1^ ■ r) (In 2 • f)] . (28) 

The interaction strength g = (47rft/io7jv) / ^ depends on the intra-distance R between two spins. The unit vector f 
connects two spins and characterizes the alignment of spin pair. If we do not apply an additional magnetic field, the 
system energy is only dependent on the value of g, thus it is not possible to obtain information about the direction 
f. On the other hand, when the interaction strength is not very strong, the amplitude of driving on NV spin (which 
should matche the target system transition frequency) will not be strong enough to suppress the effect of the 13 C 
spin bath in diamond. Therefore, we propose to apply a strong magnetic field ^nB ^> 9- Under such a condition, the 
flip-flip process is unlikely and the system Hamiltonian can be rewritten as 

h s = uj n (r Nl + r N2 ) + g 




TZ TZ ^ ( TX TX | jy jy \ 

1 N 1 1 N 2 - 9 V A^i A^2 + I N 1 I N 2 ) 



(29) 
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FIG. 9. (Color online) Eigenstates of two interacting spins under a strong magnetic field. We use NV spin to detect the transition 
frequencies of Qi and 0.2, the difference between which provides the information about the spin-spin coupling strength and the 
relative direction of the spin pair with respect to the applied magnetic held. 



where the spin operators I x ,I y ,I z are defined in the quantization axis induced by the applied magnetic field. The 
eigenstates can be written as 



\Eq) = |tt) , E = lo n + I [l - 3(f • b{6, <P)) 2 



\Ei) = \l\{\n) + \tt)),E 1 



l-Z{r-b{9,<p)f 



\E2) = \/^(\n)-\m,E 2 =0, 



\E 3 ) = \U),E 3 = -u N - 



(30) 

(31) 

(32) 
(33) 



The transition from the singlet eigen state to the states \Eq) and \Es) is determined by the difference of the 
coupling operators of two spins to the NV center, namely A\ — A^. Therefore, if we tune the Rabi frequency of the 
NV spin around the Zeeman energy cjtv, the dominant flip-flop process from the transition from \E\) to \Eq) and \Es) 
with two corresponding resonant frequencies, see Figj9j 



0,2 ~ 



3fl 
4 

4 



l-3(r-b(e,<j>)) 2 
l-3(r-b(6,<p)) 2 



(34) 
(35) 



The difference between these two resonant frequencies is 



35 



l-3(r-b(6,<f>)) 2 \, 



(36) 
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thus gives us information about the coupling strength g and the alignment vector f. To determine these exact values, 
we propose to apply magnetic fields in nine directions, and get the resonant frequencies respectively as follows 
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(37) 
(38) 
(39) 
(40) 
(41) 
(42) 
(43) 
(44) 
(45) 



After some calculations, we can obtain 
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FIG. 10. (Color online) Resonance of a single 1 H20 molecule with a NV center under an additional magnetic field 7ifj£> = 
500kHz in different directions: x (a), y (b), z (c), x + y (d), x — y (e), x + z (f), x — z (g), y + z (h), y — z (i). The signal S is 
measured at time t. In each plot, A denotes the distance between two resonant frequencies around the Larmor frequency as in 
Eq. (34-36). The molecule is 5nm from the NV center, and the vector that connects two hydrogen nuclei is described by 

(r, 0o, 0o ) = (0.1515nm, 118.2°, 288.85°). From these results, we can infer that the distance is r — 0.1518nm, and the alignment 
direction is (0, 0) = (118.29°, 288.82°) following Eq.(50-52). 
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From Eq. (37-48), we obtain the value of coupling strength g as 
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(46) 
(47) 
(48) 

(49) 

(50) 
(51) 
(52) 



Finally, with the obtained values of g 2 



' xi ' yi 



we compare Eq. (40-48,50-52) and know the relative signs between 



c ,r y ,r z and thereby we get the align direction r of the spin pair. 

In the main text, we have applied our ideas to measure the distance between two hydrogen nuclei in a water molecule. 



In Fig 10 , we plot the resonance frequencies with a magnetic field in nine different directions. The strength of the 
magnetic field is 7ifj^ = 500kHz. The amplitude of driving on NV spin is thus strong enough to suppress the effect 
of the 13 C spin bath, see Figj5] (b). The 1 H2 molecule is 5nm from the NV center, and the vector that connects 
two hydrogen atoms is described by (r, #o,</>o) = (0.1515nm, 118.2°, 288.85°). With the resonant frequencies as in 



Fig 10 , we calculate based on Eq. (49-52) and obtain that g = 34.684kHz, and thereby the intra- molecular distance is 



d = 0.1518nm, and the alignment direction is (0,0) = (118.29°, 288.82°). 



Measure distance between two organic spin labels 

Spin labels are organic molecules with a stable unpaired electron. They can be attached to the protein (covalent 
or as a ligand) via a functional group. Electron spin resonance (ESR) based on spin labels has been used as a 
spectroscopic ruler to determine protein structure and monitor macromolecular assembly processes. It is however 
hard to go beyond a distance of 5nm, due to e.g. (in) homogeneous line broadening. Here, we propose to use the 
NV center as a probe and combine with the application of spin labels, which opens possibilities of measuring a intra- 
or inter- molecular distance larger than 5nm at a single molecule level [5]. We consider the widely used nitroxide 
spin labels, in which the electron spin interacts with the nitrogen nuclear spin. The coupling is relatively strong as 
compared with the interaction between two distant spin labels. To suppress the effect of the nitrogen nuclear spin 
and enhance the signal, one can use the NV spin to first polarize the electron spin and then use the polarized electron 
spin to prepare the nuclear spin into the state \mj = 0). The coupling with the other nuclear spins is weak and can 
be suppressed by continuous driving the electron spin. If the driving amplitude is much stronger than the hyperfine 
coupling, the system Hamiltonian can be written as 



H NV _ S = tt S x + SI (S? + Sf)+g 1 - 3(f • b{6, cj>)) 



S z l S z 2 --{S x l S x 2 +SlS v 2 ) 



S Z (A 1 S Z 1 +A 2 SI), (53) 



where the quantization axis of the spin labels are induced by the applied magnetic field. In the case that ^1 ^> g, the 
system part Hamiltonian can be approximated as 



H s = Sl(Sf + Sf)-^ l-3(f-b(0,4>)) 2 SfSI 



l-3(f-6(0,0)) 2 (\0 X )(1 X \ + 11.) (0,1). 



(54) 




FIG. 11. (Color online) Resonance of a pair of nitroxide spin labels with a distance of d. The magnetic field along the vector 
which connects two spin labels. Continuous driving field with Rabi frequency of 20MHz is applied on the spin labels, (a) 
The distance is d = 5nm, the coupling strength is g = 416.1kHz, which corresponds to the difference between two resonant 
frequencies Ai = |# = 624.1kHz. (b) (a) The distance is d = 8nm, the coupling strength is g = 101.6kHz, which corresponds 
to the difference between two resonant frequencies Ai = \g — 152.4kHz. 



The eigenstates are similar to Figj9]as follows 
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The transition from the state \Ei) to the states \Eq) and \E 3 ) with two corresponding resonant frequencies 
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(55) 
(56) 

(57) 
(58) 

(59) 
(60) 

(61) 



When the distance between two spin labels is large, such that the coupling operators A\ and A2 are sufficiently 
inhomogeneous. The transition from the singlet state \E 2 ) to the states \Eq) and ^3) may also be observed with two 
corresponding resonant frequency 



^3 = SI- 



n 4 = n 
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and the difference between these two resonant frequencies is 



Ai = 



4 L 



(62) 
(63) 

(64) 



Monitor the charge recombination of radical pair 



Many chemical reactions involves radical pair intermediate, which consists of two unpaired electrons [6]. The 
radicals are in a charge separated state and interact with each other via exchange and dipole interactions. For the 
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FIG. 12. (Color online) (a) Resonance of a model radical pair reaction with a distance between two radicals of d = 2nm, the 
recombination rate is k = 1/xs -1 . The magnetic field is assumed to be along the vector which connects two radicals. The 
results are similar for other magnetic field directions. Continuous driving field with Rabi frequency of 100MHz is applied on the 
radical pair, (b) The Rabi frequency of NV spin is set at the resonant frequency. As the radical pair recombines, the effective 
flip-flop rate decreases and the slope of the signal S becomes smooth. 



distance larger than lnm, the exchange interaction is negligible, and the main contribution comes from the dipole 
interaction. Here, we consider a simple model radical pair reaction, namely the radical pair is created in the singlet 
state and recombine at the same rate k for both singlet and triplet states, the dynamics of which can be described by 
the following master equation 

j t P = -i[H, p]-^ (l ] s L sP + pL ] s L s - 2L s pL^j - X - (l ] t L tP + plS T L T - 2L T pL^ (65) 

where H is the system Hamiltonian, where the interaction between two radicals and the coupling between NV spin 
and radicals exists only when it is in the charge separated state. The Lindblad operators L^, Lt describe the 
recombination of the singlet and triplet radical pair into the product state, which are written as 

Ls = k 1 / 2 (Q S ®\P)(S\) (66) 
L T = k 1 / 2 (Q T ®\P)(S\) (67) 

with Qs and Qt are the projectors into the singlet and triplet subspace, \S) and \P) represent the charge separate state 
and product state respectively. The above master equation is equivalent to the conventional Haberkorn approach [7]. 
We apply the same idea as in the model of spin labels by applying an additional magnetic field and also continuously 
driving the radical spins to suppress the effect of surrounding nuclei. We assume the radical pair is created in the 
singlet state, so the resonant frequency near the driving Rabi frequency Q is 

(68) 

see the above section of spin labels for more details. The charge recombination leads to the decay of the effective 
flip-flop rate. If the recombination rate is comparable or smaller than the coupling between NV spin and the radicals, 
it is possible to observe the resonance frequency, see an example in Figfl2|(a). Therefore, if we tune the Rabi frequency 
of NV spin at the resonant frequency, we can monitor the recombination of the radical pairs by observing the decay 
of the flip-flop rate. In other words, it can serve as an evidence for the charge recombination. We remark that this 
is very interesting as it provides possibility to see chemical reaction process at a single molecule level and may give 
insights into how radicals recombine into the product states. 



= ft - 
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